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Abstract 

We construct numerically new axially symmetric solutions of SU(2) Yang-Mills-Higgs theory in (3 + 1) 
anti-de Sitter spacetime. Two types of finite energy, regular configurations are considered: multimonopole 
solutions with magnetic charge n > 1 and monopole-antimonopole pairs with zero net magnetic charge. A 
somewhat detailed analysis of the boundary conditions for axially symmetric solutions is presented. The 
properties of these solutions are investigated, with a view to compare with those on a flat spacetime back- 
ground. The basic properties of the gravitating generalizations of these configurations are also discussed. 

1 Introduction 

The study of topologically stable monopole 0131310] solutions to the SU (2) Yang-Mills-Higgs (YMH) equations 
with adjoint representation Higgs field, in flat spacetime, is a subject of long standing interest. Two types of 
solutions to this system are known: multimonopoles (MM) and monopole-antimonopole (MA) chains. 

The first type comprise topologically stable multimonopole solutions with topological charge n, which for n > 
\ 2 cannot be spherically symmetric 5 and possess at most axial symmetry. For vanishing Higgs self interaction 
potential, these are solutions of first order Bogomol'nyi equations and the axially symmetric multimonopole 
solutions are minimal energy and topologically stable. In the presence of a Higgs potential, i.e. when the Higgs 
coupling constant A > 0, like monopoles repel [SIE] so we would expect that the solutions of these are saddle 
points of the energy functional. Ignoringing this instability, we will include multimonopoles for YMH systems 
with A > in our definition for a MM solution. 

Recently, further to these MM solutions, new types of configurations have been considered, which represent 
composite states of monopoles and antimonopoles in (topologically unstable) equilibrium. The existence of 
such solutions was first proven by Taubes for a model featuring no Higgs potential 8 . Notwithstanding the 
absence of the Higgs potential, such solutions are not self-dual and their energies exceed the Bogomol'nyi 
bound. Kleihaus and Kunz [5| (see also JU|) constructed numerically an axially symmetric solution with zero 
net magnetic charge, consisting of a MA pair for the YMH system with and without a Higgs potential. More 
complicated configurations consisting of MA chains and vortex rings were constructed subsequently [111 IT2*] . 

The main difference between the MM and MA solutions in flat space is characterised by their distinct 
boundary conditions at infinity. They have different boundary conditions also at the origin. 

When gravity is switched on, or equivalently when the YMH system is considered on a curved background, 
the arguments leading to the saturation of a topological lower bound and hence to the first order Bogomol'nyi 
equations disappear and as a result only solutions to the second order Euler-Lagrange equations can be sought. 

Naturally, the boundary conditions employed play a crucial role in the properties of these solutions. Even in 
flat spacetime, the main difference between the topologically stable axially symmetric MM solutions on the one 
hand, and the MA type on the other, is a result of the imposition of different boundary conditions. Otherwise, 
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both types of solutions obey the same axially symmetric YMH equations . However, to date, all regular axially 
symmetric SU(2) YMH configurations of the types mentioned above have been studied on a(n asymptotically- 
) flat spacetime. It is therefore important to analyse the effects on these solutions, of different asymptotic 
structures of the spacetime. 

To this end, we are motivated to study the properties of these configurations for spacetimes with a cosmolog- 
ical constant A, in particular. In what follows we consider the case of a negative A, correspoding to an anti-de 
Sitter (AdS) geometry. This maximally symmetric spacetime has recently enjoyed a considerable amount of 
interest, motivated mainly by the proposed correspondence between physical effects associated with gravitating 
fields propagating in AdS spacetime and those of a conformal field theory (CFT) on the boundary of AdS 
spacetime ^3E]- This adds further justification to the study of classical solutions of various field theories in 
AdS, the YMH case at hand being an obvious example. 

Prominent amongst the remarkable features displayed by solutions of the gravitating Yang-Mills (YM) 
system with A < is the existence of stable particle- like and black hole solutions ^3 EH- This results from a 
particular behaviour of the YM field asymptotically, which feature is present also for the case of a fixed AdS 
background [171 118j 2 . One aim of this work is to inquire to what extent the features of YM systems with 
A < 0, persist also for YMH systems? 

Here we present numerical arguments for the existence of both MM and MA configurations, both in a fixed 
AdS spacetime background, and in the fully gravitating case. While the static solutions of the YM system are 
drastically different from those of the EYM ones, in the case of the YMH system the solutions with A < 
(with or without gravity) are found to be very similar to those of the corresponding Minkowski spacetime 
configurations. (The situation is the same also for asymptotically flat EYMH.) 

As found in |20j for A = 0, when gravity is coupled to YMH theory, a branch of gravitating solutions 
emerge from the flat spacetime configurations. This branch extends up to some maximal value of the gravi- 
tational strength. The YMH system with a negative cosmological constant has recently been studied 
for the spherically symmetric, unit magnetic charge gravitating monopole. It was found 2l : that unlike in flat 
spacetime, in a fixed AdS background there are no analytic solutions that might be used as a guide. Moreover, 
when a cosmological constant is included (no matter how small this constant is) no solution close to the flat 
space BPS configuration can be found. Numerically constructed monopole solutions have been exhibited in |22| . 
where the effects of gravity are also included. A distinctive feature of AdS solutions concerns the asymptotic 
behaviour of the fields. When A < 0, the Higgs field approaches its vacuum expectation value faster than in the 
flat space case. The radius of the monopole core decreases, as the magnetic field concentrates near the origin. 

We aim to extend this analysis for axially symmetric solutions, by studying the basic properties for both 
types of configurations, MM solutions with nonvanishing magnetic charge and MA pairs with a zero net magnetic 
charge. As expected, we find that the basic properties of the gravitating AdS configurations configurations are 
similar to the asymptotically flat spacetime counterparts. 

The paper is structured as follows: in the next Section we present the model and subject the system to axial 
symmetry. The imposition of axial symmetry is carried out in a uniform manner to include the two distinct 
boundary conditions pertaining to multimonopole and monopole-antimonopole solutions. Numerical solutions 
in a fixed AdS background are constructed in Section 3 for the two different sets of boundary conditions. Section 
4 contains a discussion of the main properties of the gravitating counterparts of these solutions. We conclude 
with Section 5 where the results are summarised. 

1 In the absence of a Higgs potential the MM solutions obey the first order Bogomol'nyi equations, while the MA solutions obey 
the second order Euler— Lagrange equations. 

2 In Ref. 1171 . a solution to the SU(2) YM system in a fixed dS background is found in closed form. Passing from dS to AdS by 
an appropriate change of sign in this solution, one discovers the same interesting feature mentioned above, which was discovered in 
115111b! for the gravitating case. 
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2 The model 



2.1 Action principle and field equations 



We consider the action for an 5*0(3) YM field A a ^ coupled to a triplet Higgs field $ a with the usual potential 
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hK» F * v + \o^ a D^ a + V(*) 



(1) 



where the field strength tensor and the covariant derivative are defined as 

F% = d^A a w - d„A« + e ahc Al Al , D^ a = d^ a + e abc A^ c , 

and g = det t^,„, g^ v being the metric of a fixed anti de Sitter (AdS) background. Subjecting the action Q to 
the variational principle results in the Euler- Lagrange equations, and varying the Lagrangian C with respect 
to the metric g^ yields the energy-momentum tensor 



T»» = F a m F^ - -g^F^g^ 6 + ]-D^ a D^ a - -g^(D a ^ a D a ^ + V(^)g^). 



(2) 



We are interested in static, purely magnetic (A t = 0), axially symmetric finite energy solutions of the Euler 
Lagrange equations. The energy density of a solution of the YMH equations is given by the tt— component of 
the energy-momentum tensor; integration over all space yields the total mass/energy 



M 
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For the mass integral to converge, each term in the integrand of (0 must vanish at large r. This will give 
the required asymptotic behaviour of the gauge and Higgs field. The additional requirements to have a finite, 
locally integrable energy density impose boundary conditions at the origin and on the z— axis. 

The static energy functional arising from is bounded from below, in flat spacetime, by the topological 
charge 

<Pxe ijk F? j D k $ a = ' 
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(4) 



which is the magnetic monopole charge, whose definition is valid also in curved spacetime. 



2.2 Imposition of axial symmetry 

We start by stating the line element for the fixed AdS background 

ds 2 = g^ u dx»dx v = dr . + r 2 {d6 2 + r 2 sin 2 6d(p 2 ) - ( 1 - -rA dt 2 , (5) 
1 - f r 2 V 3 / 

where r, 6, ip are the usual spherical coordinates. The YMH system is defined by Q), with a metric tensor g^ 
given by (JHJ). 

We proceed in this section to present the usual axially symmetric Ansatz of Rebbi and Rossi [21] > adapted 
in particular to solutions with boundary conditions supporting monopole-antimonopoles. Particular attention 
is given to finding the relevant asymptotic values. Note also that this analysis is valid in the limit A = 0. 

Expressing the SU{2) gauge connection A^ and the algebra valued Higgs field $ in terms of components 
labeled by the algebra indices 

A ll =A a J± , $ = $ Q ^ , a = a,3 = a,z ; a=l,2 = x,y, 



3 



and splitting up the spacelike coordinates x p — (2^,2:3) = (2^,2:3), with Xi — {x\,x%) = (x,y), (where r 
y/ p 2 + z 2 , 8 = arctan | ) the axially symmetric Rebbi-Rossi Ansatz can be expressed as follows 

X ' ^ (ex) in a 



A? = 


-a p 


Xi {en 


A? = 




?) 






At = 


-a z 


{en) a 



(ex)i (6) 

, 4 = 

for the gauge connection, and 

$ Q = r)(j) 1 n a , $ 3 = r 1 (f> 2 (7) 
for the Higgs field, 77 being the VEV of the Higgs field, with inverse dimension of a length. The six functions in 
©-0 depend on p — y/\x a \ 2 and 23 = z, the unit vector in the azimuthal plane is 

n a = (cos ncj), sin ncj)) , 

and eab in l|35fl is the antisymmetric Levi-Civita symbol. This ansatz contains also an integer n, which the 
winding number. 

Subject to the Ansatz © and J7J), the static energy density functional JTJ reduces to a 2 dimensional system 
in (p, z) space which exhibits a residual U(l) gauge invariance. The two components of the gauge connection 
are (a p , a z ) appearing in the Ansatz |J?jJ|, and the gauge arbitrariness can be removed by the gauge condition 

dp a p + d z a z = . (8) 

Since it is more convenient to work in spherical coordinates, we replace the two functions (a p ,a z ) by (a r ,ag) 
according to 

. cos 8 sin 8 

a„ = a r sm 8 + ag , a z — a r cos 8 — ag . 

r " r 

With this replacement, and denoting \ A = (x 1 ; X 2 ); one can write 

C = V^{g rr g ee f?g+9 vv (9 rr \VrX A \ 2 +9 ee \VgX A \ 2 ) 



n 2 



g rr \V r(p A \ 2 +g 0s \V^ A f+g^(s AB( t> A x B ) 2 ] + V(cp A ) , (9) 



which is manifestly a (7(1) gauged Higgs model, with U(l) connection (a r ,ag) and Higgs field doublets x A = 
(x 1 ^ 2 ), 4> A = [ < t )1 , c t )2 )i an d with the correspoding U(l) curvature and covariant derivatives in © defined by 

Vp X A = dp X A + apS AB x B 
V^ A = 8p^ A + a tl e AB (j ) B 

having denoted p = r,6, not to be confused with the label p, in Q. 

We now consider the asymptotic conditions in the r» 1 region that must be satisfied for finite energy, and 
the conditions of analyticity at r — and on the z-axis. 

It is immediately obvious from the Ansatz ©, in this gauge, that when p — ► 

limV-O , limv 2 = -n, (10) 

p^O p— >o 

which will be used presently to find the asymptotic values of the fields in the r>l region. In the same (p — > 0) 
limit, differentiability sets the following condition on the Higgs field function cj) 1 in Q 

lim (j) 1 oc p n , i.e. 1 | r=o = O and cj) 1 L_ „• = . (11) 

but no condition on the other Higgs field function <fi 2 . Thus the zeros of the Higgs field, i.e. when both \<p 1 \ 2 = 
and \(f> 2 \ 2 = will not necessarily occur at the origin in general. 
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2.3 Asymptotics 

In addition to the particular conditions i|l(J|) and 1111) to be satisfied on the z-axis, the comprehensive set of 
asymptotic values will now be stated, both in the r» 1 and r <C 1 regions, in that order (here also the analysis 
remains valid in the flat space limit). 

Asymptotics in the r» 1 region 

In the r > 1 asymptotic region on the other hand, the effect of a symmetry breaking Higgs potential, 
whether explicitly included in the action or not, results in the following finite energy condition on the Higgs 
field 

lim (f) 1 —smm8 , lim tfi 2 = cos m9 , (12) 

r — >oo r — >oo 

m being an integer. For m = 1, (|12J) are just the usual boundary conditions applying to (multi-) monopole 
solutions, and with n = 1 in l|7|l . they are the boundary conditions for the spherically symmetric unit monopole. 
For m > 2, the solutions describe monopole-antimonopole chains, with magnetic charges equal to the winding 
of the asymptotic Higgs field 1121) . 

Q m = 2im[l - (-l) m ] . (13) 
Finiteness of the energy requires that both terms \T> r (f> A \ 2 and \T>g(j) A \ 2 vanish at infinity, yielding 

lim a r — , lim a§ — — m , (14) 

r — >oo r — >oo 

while that for the last term in JHJ describing the interaction between \ A an d 4> A yields 

(e0) AX A = O (15) 

in this limit, where we have expressed (|12fl as 

lim 4> A = t A , (16) 

r— »oc 

in terms of the unit 2-vector t A = (sinmfl, cos m9). The 9 derivative of l|16fl gives the condition 

(et) A dex A = mt AX A , (17) 

which we will find convenient to use presently. 

The corresponding condition for the vanishing of \T> r x A \ 2 in © at infinity implies that \ A ~ * const, in 
this limit, but the actual asymptotic values of \ A cannot be determined without analysing the corresponding 
Euler-Lagrange equations, since finiteness of energy does not require that the term \T>ex A \ 2 m © vanish at 
infinity. To leading order in this limit, this equation is 

sineVg(^—Vg X \ =0. 

\ sin J 

Contracting this equation with the unit vector t A and using l|15|) and (|17fl . it reduces to 

-|m^f = o, (is) 

do sin 9 

which is immediately integrated to give 

t A dex A = cisin0, (19) 
ci being an integration constant. i|19|) an d (|15fl together now result in 

— (t A X A ) + Ci sin# = 0, 
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yeilding finally 



X A = ( c 2 + ci cos 6) t 



(20) 



c 2 being the second integration constant. 

It is easy to check that the asymptotic fields (|20|) as they stand, result in the cancellation of the singularities 
due to sin \j e in @. To fix the integration constants c\ and c 2 in 1(20(1 . additional constraints are needed. To 
this end we recall the analyticity conditions l|10(l on the z-axis. Substituting (|20|l in (|10Jl . we find that the first 
member of the latter is identically satisfied while the second yieds 

x\=a = (c 2 + Ci) = -n 
x\=, = (-l) m (c 2 -c 1 ) = -n ) 

yielding the values of c\ and c 2 to be used in (|2U)l . 

ci = , C2 = — n, to even (21) 
c\ = —n , c 2 = , to odd . (22) 

Asymptotics in the r ^ 1 region 

The r <§C 1 region is a special case of the p = line, namely the z-axis. Starting from the constraints of 
analyticity on the z-axis, we have already i|l(J|) - (|ll|) . We restate these asymptotic behavioures in the p — > 
limit, along with the other such conditions implied by the Ansatz ||§J, although we will be concerned almost 
exclusively with the actual asymptotic values. For the Higgs field funtions cf> , A = 1,2, these are 

cj> 1 = a l {z) P n , c/> 2 = a 2 (z), (23) 

and for the gauge field functions (a p , a z ) and X A A = 1, 2, 

a„ = 6i( 2 K +1 , a z = 6 2 (z)p" , x 1 =c 1 (z)p"+ 2 , x 2 = _„ + c 2 {z)p 2 . (24) 

It is clear from (I24|l that the asymptotic values of these functions are simply stated as 

a p | r=0 = , a z | r=0 = , x 1 |r=o = , x 2 | r =o = -n . (25) 

To state the asymptotic values of the functions (cj) 1 ,^ 2 ) on the other hand is much harder. It is well known that 
for axially symmetric multimonopoles |25| . with to = 1, these are 

^|r=0 = , tf\r=0 = 0, (26) 

while for monopole-antimonopole-chain solutions with m > 2, the situation is more complex. Clearly, for 
monopole-antimonopole-chains on the z-axis, the function a 2 (z) in l|23() must have zeros away from r — 0, 
precluding the second member of H26(l . To date, where these zeros occur has not been determined analytically, 
but are found through the numerical process ^2 1201 • 

For monopole-antimonopoles (|26Jl is inadequate, so from 1)23(1 we infer the general, weaker, boundary condi- 
tions 

^\ r =o = , c^ 2 | P=0 = 0. (27) 

We note that both 1(26(1 and 127(1 are consistent with the conclusions of |2E], namely that theere is only one zero 
of the Higgs field for an axially symmetric MM solution, and that if there are distinct zeros on the z-axis, either 
axially symmetry is violated or the charges corresponding to the distinct zeros have differening signs. 

The parametrisation of Kleihaus and Kunz 

This particular parametrisation of the Rebbi and Rossi Ansatz used by these authors, see e.g. j2Zj) is very 
convenient and is employed in most works on axially symmetric YMH systems. Since the numerical integrations 
will be carried out for the variables (r, 9) and not (p,z), it is useful to employ a parametrisation for which 
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the Ansatz ©-0 agrees with the unit magnetic charge monopole for n = 1. The functions (Hi, H 2 , H 3 , H4) 
parametrising the gauge field are related to the functions in the axially symmetric Ansatz © as follows 

Hi = r(a p sin 9 + a z cos 9) 
1 — Hi = r(a p cos 9 — a z sin 9) 

(28) 

n sin 9H 3 = x 1 sm # + ( rl + X 2 ) cos # 
nsin#(l — H4,) = x 1 cos^ - (n + x 2 ) sin0 , 

and the functions ($1, $2) parametrising the Higgs field Ansatz Q, 

$ x = <^ x sinff + 2 cos , $ 2 = ^ cos 9 - 2 sin 9 . (29) 

In terms of this parametrisation, the boundary conditions for the monopole-antimonopole solutions at 
r — > 00, can be read off (fTljl . |2"0" )l -l|2"2" )l . and J2HJ- For odd m these are 

COS $ COS $ 

ffi=0, ff 2 = -(m-l), #3 = -^r[cos(m- 1)9- 1], ff 4 = — — sin(m-l)0. (30) 

sm# smfci 

and for even m 

Hi = 0, H 2 = -(m-l), # 3 = ^-[cos(m-l)0-cos0], H 4 = _ sm ( m - V e (31) 

sm sm t/ 

while 

$1 = cos(?ti - 1)9, $ 2 = sin(m - 1)0, (32) 

for any value of m. 

The corresponding boundary conditions for multimonopoles are simply given by the set for odd m, 130(1 . 
with m = 1. 

The boundary conditions for the gauge field functions at r — > can be read off 1(25(1 . for both multimonopoles 
and monopole- antimonopoles. The Higgs field boundary conditions for MM's and MA's on the other hand, are 
different, being stated respectively by (|2"tJl) and i(2"7|) . Using the relations lt2"51)- lt2"9T) . (|2~5j) and l(2l)|) yield the 
required boundary values at r = for multimonopoles 

Hi=H 3 = 0, H 2 =H i = 1, $1 = $ 2 = 0. (33) 

For monopole- antimonopoles, the less stringent condition 1(27(1 must be used instead of ((26(1 . Following |20llll| 
we state the boundary values for monopole- antimonopoles as 

Hi=H 3 = 0, H 2 = Hi = l, cos0d r $ x -sin0d r $ 2 = 0, sin 6> $ x + cos0 $ 2 = , (34) 

the Higgs part of which is a weaker version of ((27(1 . 

The above parametrisation is in fact employed in the literature for the numerical construction of axially sym- 
metric MM solutions, e.g. in |22] and also for gravitating YM solutions 28 . In the corresponding construction 
of axially symmetric MA solutions however, these authors use a different parametrisation in Ref. j^j, as well as 
in |20j . where they have constructed the gravitating counterparts of the latter. However, the parametrisation 
used in the MA works |2U] is equivalent to that of |J7], which is what we have finally used in the present 
work, unifying the treatments of axially symmetric MM and MA solutions. 
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Figure 1. The mass of solutions is plotted as a function of the cosmological constant for three values of the wind- 
ing number n and Higgs self-coupling constant A = 0.5. Here and in Figures 2 and 4a the mass is given in units 
4-Krj/g. 

3 Solutions in a fixed AdS background 

We have performed our numerical computations using the parametrisation of Kleihaus and Kunz described 
above. All numerical calculations were performed with the software package CADSOL/FIDISOL, based on the 
Newton- Raphson method |29j . 

3.1 Multimonopole solutions 

We start by discussing the simplest example of axially symmetric solution, describing multimonopoles containing 
only magnetic charges with the same sign (m = 1). The suitable boundary conditions for the Higgs field at 
infinity are <&i = 1, $2 = 0, while, in the same limit the YM potentials should vanish Hi = 0. 

Given the parity reflection symmetry satisified by the ansatz H28[) - <|29l) . we ne ed to consider solutions only 
in the region < 9 < tt/2; on the z- and p-axis the functions H\,H^,^2 and the derivatives deH2,dgH4 and 
are to vanish. To fix the residual abelian gauge invariance we choose the usual gauge condition JSJ which 
reads in terms of Hi functions 

rd r Hi - d e H 2 = 0. (35) 

For these boundary conditions, the configurations with n = 1 corresponds to spherically symmetric monopoles 
with Hi = H3 = $2=0, H2 = H4 = w(r), $1 = H(r), and generalize the 't Hooft-Polyakov solution for a 
negative cosmological constant. Axially symmetric solutions are found for n — 2,3, ... . 
Dimensionless coordinates and Higgs field are obtained by rescaling 

r— >r/r), A^Ary 2 , $^77$. (36) 

Within this ansatz and gauge choice, we have solved the resulting set of six non-linear coupled partial 
differential equations, finding solutions for any considered value of the cosmological constant. 
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Figure 2. The mass per topological charge is plotted as a function of Higgs self-coupling constant A for MM solu- 
tions with A = —0.5, — 2 and n = 1, 2. 





Figure 3. The dimensionless energy density of a n = 2 MM solution in a fixed AdS background with A = — 1 is 
plotted as a function of p and z. The solution has been found for a Higgs coupling constant A = 10 and has a mass 
M = 3.998. 

Qualitatively, the behaviours of the Higgs field and Yang-Mills fields are very similar to that corresponding 
to Minkowski spacetime monopolcs. In particular, we notice a similar shape for the functions Hi and $i and 
also for the energy density. Again, the functions H 2 , H4 and $1 present a small 9— dependence. For this type 
of solutions, the nodes of the Higgs field are superposed at the origin. 

In Figure 1 the mass of the solutions in units of ^nrj/g is plotted as a function of the cosmological constant 
for various winding numbers and a fixed value of the coupling constant A = 0.5. 

The energy of the multimonopole is of the order of n times the corresponding one-monopole energy. However, 
for any value of A, the ratio M/n increasing with increasing n. Also, this quantity is greater than the mass of 



the corresponding n = 1 spherically symmetric monopoles. For example, for solutions with A = —4 and A = 0.5, 
we find M(n = 2)/M(n = 1) = 2.466, M(n = 3)/M(n = 1) = 4.373 while M(n = 4)/M(n = 1) = 6.678. 

Most of our results have been obtained in the Prasad-Sommerfeld limit, but similar conclusions hold for a 
nonvanishing Higgs self-coupling constant A. To demonstrate this dependence, we plot in Figure 2 the mass 
of solutions per topological charge as a function of A for two different values of the cosmological constant and 
winding numbers n = 1, 2. 

For all configurations, the energy density e = — T t of the solutions has a strong peak along the p axis, and 
it decreases monotonically along the symmetry axis. Individual unit charged monopole components of the MM 
do not feature distinct peaks (see Figure 3). Equal density contours reveal a torus- like shape of the solutions. 
For increasing |A|, these peaks increase in size and the energy density becomes localised in a decreasing region 
of space. 

3.2 Monopole-antimonopole solutions 

As discussed in [H] for n = 1 and A = 0, there exist also a different type of solutions of the second order 
Euler-Lagrange equations, which are not stable and represent saddle points of the energy, rather than absolute 
minima resulting from solving the Bogomol'nyi equations. Here we construct the simplest examples of such 
solutions, namely MA solutions obtained by taking to = 2 in the general Higgs field asymptotics at infinity, 
(|31|I - H32(I . which implies the asymptotic behaviour ifi = if 3 = 0, if 2 = H4 = — 1. This corresponds to a MA 
configuration, with a vanishing net magnetic charge. 

Similarly to the multimonopole case, we need to consider solutions only in the region < 9 < n/2; on the 
z- and p-axis the functions ifi, if^, $2 and the derivatives dgH2, doH^ and dg&i are to vanish. The same gauge 
condition, l|35|l . fixes the residual Abelian gauge invariance. 

A calculation similar to that done in [Hj shows that this configuration possesses two magnetic charges of 
opposite sign, located on the positive and negative z— axis, respectively, with nonvanishing local density of 
magnetic charge. However, the integral (@J, evaluated at infinity, for a surface enclosing both charges, yelds a 
zero net magnetic charge. Also, the expansions of the matter functions at the origin and on the z—axis presented 
in Ref. [5] are still valid, and will not be given here. 

The magnetic dipole moment C m of these solutions can be obtained from the asymptotic form of the non- 
abelian gauge field, after choosing a gauge where the Higgs field is asymptotically constant $ — > T3, which 
yields 

A k dx k = Cj—^dip. (37) 
r 2 

We have constructed MA solutions for a large range of the parameters (A, A). For vanishing A, our results 
are in very good agreement with those of UJ. The main part of the numerical analysis has been done for the 
case n = 1, the case n = 2 being studied more briefly. However, most of the qualitative properties of the n = 2 
solutions do not differ from those of the n = 1 case. 

Similarly to the case of multimonopoles, we have found that a negative cosmological constant does not 
change the qualitative properties of the solutions. As expected, a negative cosmological constant affects the 
behaviour of the fields in the asymptotic region, where the modulus of the Higgs field, for example, reaches its 
asymptotic value faster. 

In Figure 4 we present the the mass and the magnetic dipole moment of n = 1 MA solutions as a function 
of A for three values of the Higgs self-coupling constant A. 

The energy density always possesses maxima on the positive and negative z-axis at the locations of the 
monopole and antimonopole and a saddle point at the origin. An increase of |A| makes the maxima of the 
energy density higher and sharper. At the same time, the modulus of the Higgs field tends faster towards its 
vacuum expectation value. The typical energy density for nongravitating MA solutions have a similar form to 
that presented in Section 4 for gravitating configurations. 
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Figure 4. The mass and the magnetic dipole moment of n — 1 MA solutions is plotted as a function of A for three 
values of the Higgs self-coupling constant A. 

The nodes of the modulus of the Higgs field correspond to the locations of the particles. With increasing 
|A|, the distance d between the MA centers becomes smaller tending to a finite limit as A — > -co, making it 
more difficult to distinguish the monopole from antimonopole. Several values we found for solutions with A = 2 
are: d(A = 0) = 3.24, d(A = -0.5) = 1.52, d(A = -1) = 1.22, d(A = -2.5) = 0.9, the distance between the 
MA centers approaching for larger values of |A| a limiting value d ~ 0.85. The magnetic dipole moment C m 
decreases in the same limit. 

Considering solutions with a fixed negative value of A, we find that the distance between the MA centers 
decreases as A increases and converges to a finite, nonzero value. For example, for MA solutions in an AdS 
spacetime with A = —1, d decreases from d(X = 0) = 0.68 to a limiting value d ~ 0.6. This property is shared 
with the A = case [5]. 
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4 Inclusion of gravity 



We now consider the effects of these axially symmetric configurations on the AdS spacetime, by including the 
Einstein gravity term with a negative cosmological constant y/— g{— R/(16t:G) + 2A) in the action density in 
0. 



The usual rescaling l|36|) reveals in this case the existence of a new dimensionless coupling constant 



2 



AnGr/ 2 , where G is the Newton constant. The complete classification of the solutions in the space of physical 
parameters (a, A, A) is a considerable task which is beyond the scope of this paper. Instead, we analysed the 
situation for several values of A which, hopefully, reflect all the properties of the general pattern. Also, for 
simplicity we will study in this section only solutions with no Higgs potential (A = 0). 

The asymptotically anti-de Sitter (AAdS) gravitating counterparts of the axially symmetric configurations 
discussed in Section 3 are found by using a metric of the form 

d * 2 = T^T^TT + r2d92 ) + si * 2 0d ^ - /C 1 ~ (38) 
J 1 — j r 2 J 6 

with /, I and m being functions of r and 9. To obtain asymptotically AdS regular solutions with finite energy, 
the metric functions have to satisfy the boundary conditions / = m = Z = lat infinity, and d r f = d r m — d r l = 
at the origin. The boundary conditions on the z-axis are dgf — dgm — dgl = 0, which are consistent with the 
requirement of invariance under reflection in the 9 = tt/2 plane. Finally, regularity on the z-axis requires also 
in == I for 9 = 0. These metric boundary conditions are valid for both MA and MM gravitating configurations. 
The boundary conditions for the gauge potentials and Higgs field are similar to the nongravitating case. 

To solve the set of nine EYMH equations we employ the same numerical algorithm as for the YMH solutions 
in fixed AdS background. In the numerical procedure we use a suitable combination of the Einstein equations 

Rfiv - 2^" R + = 87rGT^, (39) 

such that the diferential equations for metric variables (f,l,m) are diagonal in the second derivatives with 
respect to r. 

Similar to the asymptotically flat case, the value of the mass is encoded in the asymptotics of the metric 
functions f,l,m, whose expression, valid for large r, is 



A + f 2 sin 2 9 1 mi +m 2 sin 2 9 1 h + 1 2 sin 2 9 1 
J = H 3 l-O(-r), m = H = + 0(-r), i = H 3 + 0{-r), (40) 



where /i, fa are constants to be determined numerically, and 

= = ~Y' = 17' 171,2 = ~vf' 

The generalization of Komar's formula for AdS asymptotics is not straightforward and requires the further 
subtraction of a background configuration in order to yield a finite result. 

To compute the mass of these configurations, we employ the counterterm formalism proposed by Balasub- 
ramanian and Kraus |33| . This technique was inspired by AdS/CFT correspondence and consists of adding 
suitable counterterms I ct to the total action of the theory. These counterterms are built of curvature invariants 
of a boundary metric and thus obviously do not alter the bulk equations of motion. However, they do regularise 
the boundary stress tensor and also the conserved charges. The details of the computation, which we omit here, 
are presented in |34j . where an axially symmetric line element with the same paramctrisation and the same 
asymptotics is considered. 

The expression of the mass derived in this way is 

The same form is obtained by using the Hamiltonian formalism of Henneaux and Teitelboim |32| . 
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Figure 5. The metric function / and the norm of the Higgs field $ = of typical AAdS MM solutions 

are shown as a function of r for three values of the angle 9 with n = 2, A = —3 and a sequence of values of a. 



The dimensionless mass is given by 

M /a 2 = y—M. (43) 

47T?7 

4.1 Gravitating MM 

A spherically symmetric line element is obtained for I = m, with / and m being functions of the coordinate r only. 
Some properties of the corresponding n = 1 monopole solutions are discussed in |22| by using Schwarzschild-like 
coordinates. Similar results are found also by using the above metric parametrisation. 

As found in |22j . when a is increased from zero, while A is kept fixed, a gravitating monopole branch emerges 
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Figure 6. The mass per topological charge is plotted as a function of a for n = 1, 2, 3 (multi-)monopole solu- 
tions with A = — 10 and A = 0. 

smoothly from the corresponding AdS space monopolc solutions. The total mass of gravitating solutions de- 
creases with increasing a. A similar property has been noticed |23) for asymptotically flat solutions too. 

The fundamental n — 1 monopole branch extends up to a maximal value cx max of the coupling constant 
a. This maximal value decreases with |A|; for example a max (A = 0) ~ 1.4, a max (A = —0.1) ~ 1.28, while 
a mnl (A = —3) ~ 0.80. As is known from the work of [23], Ior A = 0, and with the Higgs coupling A = 0, 
the fundamental monopole branch bends backwards at a max ^ until a critical value a cr corresponding to the 
extremal Reissner- Nordstrom (RN) solution is reached. Fro A < 0, we refer to such RN solutions as RNAdS. 
However, we have found that for large enough values of |A|, the maximal value a max and the critical value a cr 
seem to coincide, as they do also for the asymptotically flat case when A > ■ 

Along the fundamental branch, the metric function f(r) develops a minimum, which decreases asymptotically 
to zero, as a — ► a cr . The functions uj(r) and H(r) parametrising the gauge and the Higgs fields respectively, 
approache their respective RNAdS values u> = 0, H = 1 [22]. It is also noticed that the critical value of a 
as well as the corresponding value of mass decrease with A. The limiting spacetime consists of an inner part 
with r < r c and an outer part with r > r c . The exterior of the critial solution corresponds to the exterior of 
a RNAdS black hole solution with a degenerate horizon at r = r c and unit magnetic charge. A discussion of 
these solutions and typical profiles are presented in [ST1I22] . 

Here we have constructed axially symmetric multimonopoles for n > 1, the metric functionsin this case 
displaying an angular dependence. Most of the properties of these solutions are similar to those corresponding 
to asymptotically flat space. In particular, for a given topological charge n, we find a branch of globally 
regular asymptotically AdS multimonopoles, emerging smoothly from the corresponding solutions in fixed AdS 
background. 

In Figure 5 we show the metric function / and the norm of the Higgs field |$| = + °f gravitating 
n = 2 multimonopolc solutions with A = —3, for several values of a. 

The fundamental branch extends up to some maximal value of a, whose precise value depends on A. Along 
this fundamental branch, the mass of the solutions decreases monotonically. For the studied case A = 0, the 
value of a max decreases with increasing A. The value at the origin of the metric function / decreases with 
increasing a, corresponding as a — > a cr to the degenerate horizon of a RNAdS solution. At the same time 
the Higgs field function <\>\ approaches the unit value. The exterior of the critical solution corresponds to the 
exterior of a degenerate horizon RNAdS black hole with magnetic charge n (note that for the same value of a 
and rii > ni, the mass of this RNAdS black hole satisfies the relation Min-i) > M{n{f). 
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Figure 7. The dimensionless energy density of a n — 1 gravitating MA solution with A = —0.04, a = 0.2, is plotted as 
a function of p and z. 

The mass per unit charge of the multimonopole solutions decreases with increasing a. As discussed in \HJ\, the 
mass per unit charge n of the A = multimonopoles is smaller than the mass of n = 1 monopoles, i.e. that like 
monopoles attract. We have found however that this is not a generic property of AAdS solutions, and that for 
large enough values of A the monopoles of like charge in a gravitating Higgs model on AdS space repel for all 
values of a (a typical situation is presented in Fig. 6). 

4.2 Gravitating MA 

For completeness, we present here some properties of the gravitating MA solution with A < 0. The field 
equations have been solved in this case by using the same metric ansatz (13811 . and taking again m = 2 in the 
matter field asymptotics at infinity, H31fl - I|32l) . for n = 1, A = and several values of A. 

For MA boundary conditions, we find that, similar to the MM case, a branch of gravitating MA solutions 
emerges smoothly from the MA solutions in fixed AdS background. The basic properties of these solutions are 
similar to those of the configurations in a fixed AdS background. For example, the modulus of the Higgs field 
always possesses two zeros on the z— axis, corresponding to the location of the monopole and the antimonopole. 

With increasing a, the monopole and antimonopole move closer to the origin and the mass p of the solutions 
decreases. In Figure 7 we plot the energy density e = — T\ of a typical MA solution as a function of the 
coordinates p, z, for A = —0.04, a = 0.2. 

This fundamental branch of solutions ends at a critical value a cr , when gravity becomes too strong for 
solutions to persist. As expected, this value decreases with |A| (for example a cr (A — 0) ~ 0.67 (201, cn cr (A = 
-0.01) ~ 0.64 while for A = -1 we find a cr ~ 0.35). 

In asymptotically flat space, it has been found that at a cr a second branch of MA solutions emerges, extending 
back to a — [201 . Along this upper branch, the distance d between the MA centers shrinks to zero size in 
the limit a — > 0. After a suitable rescaling, as a — * the configuration approaches the spherically symmetric 
Bartnik-McKinnon (BK) solution, with Hi = H% = 0, H2 = H4 — u>(r), with u>(r) — > — 1 asymptotically (thus 
a — > means here 77—5-0 with fixed G). In AdS spacetime however, the situation is to be more complicated 
since the BK configuration is now replaced by a continuum of solutions. Here the picture depends crucially on 
the value of A (in particular, solutions with u> — > —1 cease to exist for a large enough |A|), and the meaning of 
the limiting solutions is less clear. We hope to come back to this point in future. 
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5 Conclusions 



We have studied, analytically and numerically, the basic properties of two types of solutions to the SU(2) 
Yang-Mills-Higgs model in the presence of a negative cosmological constant. The first type of solutions studied 
are the axially symmetric multimonopoles (MM) with magnetic charge equal to the azimuthal winding number 
n. The second type are the axially symmetric monopole-antimonopole (MA) pairs with vanishing magnetic 
charge and azimuthal widing number n, and nonvanishing magnetic dipole moment. These are the simplest 
examples of MA chains and vortex lines, and are not topologically stable even in flat spacetime, unlike the MM 
solutions. Both MM and MA solutions obey the same Euler-Lagrange equations of this system, subjected to 
axial symmetry. The drastic difference in the said properties arises from the imposition of different boundary 
conditions at infinity, in each case, and while the MM solutions are labeled by an azimuthal winding number 
n, the MA solutions are further labeled by a number m multiplying the boundary value of the polar angle, i.e. 
by the pair of integers (n, to). We have emphasised this point and have presented a unified treatment of both 
types of solutions. 

The spherically symmetric analogues of the MM solutions studied here were discussed in [22] , where unit 
magnetically charged monopoles on both fixed AdS backgrounds and gravitating ones were considered, but not 
MA type solutions which are at most axially symmetric. 

Qualitatively, the behaviour of the AdS solutions we constructed were found to be very similar to that of 
corresponding Minkowski spacetime configurations. Thus, it seems that when studying gauge field systems 
containing a scalar field, the Higgs field in the case at hand, the solutions exhibit a generic behaviour for A < 0. 
A similar behaviour has been noticed for sphalerons |30| . 

Let us list some points of contrast and similarity of MM and MA solutions, between the A = and A < 
cases: 

• With |A| > 0, the MM solutions satisfy the second order Euler-Lagrange equations and not the first order 
Bogomol'nyi equations, even when the Higgs potential is absent, i.e. when A = in The immediate 
consequence of this is that in this limit (A = 0), like charged monopoles are not non-interacting. It turns 
out that the mass of one n = 2 MM is larger than that of two n = 1 MM's with their centres infinitely far 
apart. Thus, it looks as if the n — 2 axially symmetric MM is unstable since like monopoles repel each 
other. This feature holds for all n, and is in contrast to the MM solutions in flat spacetime, which are 
stable for any n since they are all self-dual solutions. 

• When |A| > 0, it turns out that like monopoles remain mutually repulsive even for positive Higgs coupling 
constant A > model. This property, which is not surprising, is similar to the flat spacetime, A = case. 
Like in the case of flat background moreover, when A ^> 1, the mass of the MM solution tends to a finite 
limiting value which we have not estimated here. This may be an interesting detail to return to in 
future. 

• It is known from 9 that in the flat space model the distance on the z-axis between the monopole and 
the antimonopoled ecreases as the Higgs coupling constant A increases, converging to a finite limit. In the 
present model with |A| ^ 0, this property is similar for the MA solutions. For a fixed negative value of 
the cosmological constant, this distance decreases as A increases and converges to a finite value. This is 
qualitatively the same as for the flat space model, but quantitatively this distance decreases faster with 
increasing A. 

• For a fixed value of the Higgs coupling constant A, increasing |A| ^ of the negative cosmological constant 
results in the decrease of the distance between the monopole and the antimonopole of the MA solution, 
converging to a finite value. 

• Another one of the properties of in Minkowski space that persists in the AdS case, is that found by 
Houston and O'Raifeartaigh [21>|: Any regular axially symmetric magnetic charge distribution can be 
located only at isolated points situated on the axis of symmetry, with equal and opposite values of the 
charge at alternate points. In particular, if only one sign of the charge is allowed, all the charge must be 
concentrated at a single point. The results of the present paper support an extension of this property to 
the AdS case. 
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As expected, the inclusion of gravity does not change the general picture, since for A < 0, the gravitational 
field does not affect the behaviour of the YMH system at infinity. As it happens in asymptotically flat case, a 
critical value for the Newton constant exists above which no regular solution can be found. However, as in the 
spherically symmetric case, this critical value is smaller than the corresponding value of the asymptotically flat 
one. 

For asymptotically flat solutions, the inclusion of gravity allows for an attractive phase of like monopoles 
not present in flat space. Here we have presented numerical arguments that for large enough values of |A|, only 
a repulsive phase exists for like monopoles. 

Hairy black hole solutions of EYMH theory with a negative cosmological constant, generalising the asymp- 
totically flat configurations discussed in |191 l2"7] can also be constructed numerically, but this was not done in 
the present work.. 

One may ask about the possible relevance of these solutions within the holographic principle and its 
AdS/CFT correspondence realisation. A scalar field has been discussed in this context by many authors. 
For example, the asymptotic behaviour of the bulk scalar field was directly related to one-point functions in the 
dual CFT The vortex solution of the Abelian Higgs model has also an interesting AdS/CFT interpretation, 
its mass density being dual to the discontinuity of the logarithmic derivative of the correlation function of the 
boundary scalar operator |36| . To date however the crucial interplay between the non Abelian gauge fields and 
a scalar multiplet (with its rich set of boundary conditions discussed in this paper) is not considered in this 
context in the literature. 

On the other hand, the EYMH bulk action (or a suitable modification supporting the same set of asymptotic 
boundary conditions) does not seem to correspond to a known supergravity truncation, and in particular we do 
not know the underlying boundary CFT. We believe that further work in this direction will be of interest. 
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